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Abstract 

We study the complex-temperature phase diagram of the square-lattice Ising model for 
nonzero external magnetic field H, i.e. for < \i < oo, where \i = e~ 2/3H . We also carry 
out a similar analysis for — oo < /i < 0. The results for the interval — 1 < /i < 1 provide 
a new way of continuously connecting the two known exact solutions of this model, viz., at 
fi = 1 (Onsager, Yang) and /i = — 1 (Lee and Yang). Our methods include calculations of 
complex-temperature zeros of the partition function and analysis of low-temperature series 
expansions. For real nonzero H, the inner branch of a limagon bounding the FM phase 
breaks and forms two complex-conjugate arcs. We study the singularities and associated 
exponents of thermodynamic functions at the endpoints of these arcs. For fi < 0, there 
are two line segments of singularities on the negative and positive u axis, and we carry out 
a similar study of the behavior at the inner endpoints of these arcs, which constitute the 
nearest singularities to the origin in this case. Finally, we also determine the exact complex- 
temperature phase diagrams at fi = — 1 on the honeycomb and triangular lattices and discuss 
the relation between these and the corresponding zero-field phase diagrams. 
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1 Introduction 



The two-dimensional Ising model serves as a prototype of a statistical mechanical system 
which undergoes a phase transition with associated spontaneous symmetry breaking and long 
range order. The free energy of the (spin 1/2) Ising model was first calculated by Onsager 
JU, and the expression for the spontaneous magnetization first derived by Yang (both 
for the square lattice). However, the model has never been solved in an arbitrary nonzero 
external magnetic field, and this has long remained an outstanding open problem. Hence, 
any additional information that one can gain about the Ising model in a magnetic field is 
of value. In elucidating the properties of the zero-field model, it has proved useful to gener- 
alize the temperature variable to complex values. There are several reasons for this. First, 
one can understand more deeply the physical behavior of various thermodynamic quanti- 
ties by seeing how they vary as analytic functions of complex temperature (CT). Second, 
one can see how the physical phases of a given model generalize to regions in appropriate 
complex-temperature variables. Third, a knowledge of the complex-temperature singularities 
of quantities which have not been calculated exactly helps in the search for exact, closed- 
form expressions for these quantities. For the (spin 1/2) zero-field Ising model on the square 
lattice the complex-temperature zeros of the partition function were first discussed in Refs. 
||. [1| and the associated phase diagram is known exactly. In the complex Boltzmann weight 
variable z = e~ 2K (see below for notation) the phase boundaries consist of two intersecting 
circles with a z — > —z symmetry; more compactly, in the variable u = z 2 which incorporates 
this symmetry, they are given by a limagon of Pascal ||. Just as this generalization to 
complex-temperature has yielded a deeper insight into the zero-field model, so also it can 
shed light on the behavior of the model for nonzero field. Accordingly, in this paper, we 
shall investigate the complex-temperature properties of the 2D Ising model for nonzero ex- 
ternal field. Our methods include calculations of complex-temperature zeros of the partition 
function and analyses of low-temperature series expansions. 

Although no one has solved the 2D Ising model in an arbitrary field, Lee and Yang did 
succeed in solving exactly for the free energy and magnetization for a particular manifold 
of values of H depending on the temperature T, given by H = i(n/2)k B T (see also 
0). Although this is not a physical set of values, owing to the imaginary value of H and 
the resultant non-hermiticity of the Hamiltonian, this model is nevertheless of considerable 
interest for the insight which it yields into the properties of the Ising model in the presence 
of a symmetry-breaking field. By an extension of our analysis to complex values of H, we 
are able to continuously connect the two known exact solutions of the 2D Ising model, at 
H = |, | and at H = i(it/2)k B T [§. 



1 



Before proceeding, we mention some related work. Rigorous results on the behavior of 
the model for nonzero external field H include the theorem that, for ferromagnetic (FM) 
spin-spin coupling, J > 0, the free energy F(T, H) is an analytic function of temperature 
0-0. In the case of antiferromagnetic (AFM) coupling, J < 0, this is not the case; there 
is a temperature Tb(H) such that as T decreases through Tb(H), the system undergoes a 
(first-order) transition from FM to AFM long-range order. Although there is no known exact 
expression for Tb(H), accurate numerical values are known (e.g., || and references therein). 
In addition to the works noted above [§, [§] , papers on complex-temperature singularities in 
the 2D square-lattice Ising model for H = include Refs. [fL^|-|fl6| . A useful connection 



between the Lee- Yang solution for f3H = in/2 and a modified zero-field Ising model with 



certain coupling(s) shifted by iir/2 was discussed in Refs. [|17|, jl8| . The complex-temperature 
phase diagram for the Lee- Yang case f3H = in/ 2 was worked out in Ref. inlcuding 
exact results for the specific heat and magnetization critical exponents and series analyses to 
determine the susceptibility critical exponents at certain complex-temperature singularities. 
Complex-temperature zeros of the partition were studied for the 3D Ising model in a nonzero 
magnetic field in Ref. p0 |. In contrast to the Ising model, a certain superexchange model 



can be solved exactly in a field |2~l l. We have studied the complex-temperature properties 



of this model for nonzero field and will report the results elsewhere. Whereas here we study 
complex-temperature zeros of the partition function for real and certain complex values of 
magnetic field, a different complexification is to study the zeros of Z in the complex magnetic 
field for physical temperature |6| . This led the famous theorem || that for ferromagnetic 
couplings the zeros in the complex e~ 2 @ H plane lie on a circle and pinch the positive real 
axis as the temperature decreases through its critical value; this work also derived integral 
relations expressing thermodynamic quantities in terms of integrals of densities of complex- 
field zeros. Similar relations for complex-temperature zeros were later discussed in Refs. 



2 General Properties 
2.1 Model 



Our notation is standard and follows that in our earlier works (e.g., ||, |ljj), so we review it 
only briefly. The (spin 1/2, isotropic, nearest-neighbor) Ising model on the square lattice is 
defined by the partition function 

Z=J2 ^ m (2.1.1) 
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with the Hamiltonian 

H=-J J2 Vn<Tn>-Hj2<7n (2.1.2) 

<nn'> n 

where a n = ±1 are the Z2 spin variables on each site n of the lattice, (3 = (/c^T) -1 , J is the 
exchange constant, < nn' > denote nearest-neighbour sites, and the units are defined such 
that the magnetic moment which would multiply the H J2 n a n is unity. (Hereafter, we shall 
use the term "Ising model" to denote this model unless otherwise indicated.) For H = 0, 
the symmetry group of the theory is Z 2 . We use the notation 

K = /3J, h = (3H (2.1.3) 

z = e~ 2K , u = z 2 = e~ AK , v = timhK (2.1.4) 

and 

fi = e~ 2h (2.1.5) 

The reduced free energy per site is / = —f3F = liniTv^oo iV" 1 In Z in the thermodynamic 
limit, where N s is the number of sites on the lattice. The zero-field susceptibility is \ = 
&M qh^ \h=Qi wnere M(H) denotes the magnetization. It is convenient to deal with the reduced 
quantity % = P~ l X- 

A useful property is that the partition function Z is a generalized polynomial (with both 
negative and positive integral powers) in u and /x. (For a lattice with odd coordination 
number, Z would be a generalized polynomial in z and fi.) On a finite lattice, for fixed /i, 
Z thus has a certain set of zeros in the u plane. Experience with the zero-field model shows 
that in the thermodynamic limit, these merge together to form curves (including possible 
line segments) across which the free energy is non-analytic. These are the only singularities 
of /, except for the trivial singularities when \K\ = 00 (which are isolated singularities and 
thus are not important for the discussion of phase boundaries.) Hence, the calculation of 
zeros of the partition function on finite lattices serves as a valuable means by which to gain 
information about the above continuous locus of points where the free energy is non-analytic 
in the thermodynamic limit .Q This continuous locus of points includes the phase boundaries 
of the complex-temperature phase diagram. It may also include certain arcs or line segments 
which protrude into and terminate in the interior of some phases and hence do not separate 
any phases. 

1 Parenthetically, we note that for anisotropic spin-spin couplings, the complex-temperature zeros would 
merge to form areas instead of curves in the thermodynamic limit fl22| . Indeed, for the (heteropolygonal) 
4 • 8 2 lattice, even if the couplings are isotropic, the zeros still form areas in this limit p3|. 



3 



2.2 Symmetries 

We record here some basic symmetries which will be used in our work. First, because Z is 
a generalized polynomial in //, in the analysis of the phase diagram, it suffices, with no loss 
of generality, to consider only the range 

-%< Im(h) < | (2.2.1) 

Second, the summand of the partition function is invariant under the transformation h — > —h, 
a n — > — a n . The sign flip h — > — h is equivalent to the inversion map 

H -> - (2.2.2) 

Hence, in considering nonzero real h, one may, with no loss of generality, restrict to h > 0. 
More generally, in considering complex h, one may, with no loss of generality, restrict to the 
unit disk 

M < 1 (2.2.3) 

in the ji plane. We shall concentrate here on the range of /i values which connect the Onsager 
and Lee- Yang solutions of the model, viz., 

- 1 < p < 1 (2.2.4) 

By the above symmetry, these values suffice to describe the entire real line 

- oo < fi < oo (2.2.5) 

Third, concerning the symmetries of the locus of points across which the free energy is 
non-analytic and the associated complex-temperature phase diagram (and, for finite lattices, 
the set of zeros of Z), these are invariant under i.e. under reflection about the 

horizontal, Re(u) axis. Fourth, if and only if ji = ±1, this locus of points is invariant under 
the inversion map 

1 

u — > — 2.2.6 

u 

This symmetry holds because of the bipartite nature of the square lattice, and this mapping 
interchanges the uniform and staggered magnetizations. Thus, if one starts at a point in the 
FM phase, the mapping takes one to a corresponding point in the AFM phase. Since, in 
general, the FM phase occupies a region in the neighborhood of the origin in the u (or z) 
plane, this shows that the AFM phase will be the outermost phase, extending to complex 
infinity in these two respective planes. This result is obvious for h = and can be seen easily 



4 



for the case /i = — 1 from the relation ]17|, |1S | connecting this theory to a zero- field Ising 



model with certain couplings shifted by i7r/2, as discussed in [[HJ. Although this inversion 
symmetry does not hold for p 7^ ±1, one can still make the following statement. Let us 
write u in polar form, u = p u e t0u . Then 

K = -~ln« = --(lnp u + id u + 2mn) (2.2.7) 

where n indexes the Riemann sheet of the logarithm and will not be important here (we take 
it to be n = 0). The usual limit of infinitely strong antiferromagnetic spin-spin exchange 
coupling K = [3 J — > — 00 corresponds to u going to infinity along the positive real axis. 
However, eq. ( |2.2.7|) shows that as \u\ — > 00 along any direction, not necessarily 9 U = 0, 
K corresponds to an infinitely strong AFM coupling, since the term — (1/A)i9 u is finite and 
becomes negligible in this limit. It follows that for arbitrary fixed p, if an AFM phase exists 
at all (it is absent, e.g., for the isotropic Ising model on the triangular and Kagome lattices), 
then it extends outward in all directions to complex infinity in the u or z plane. This general 
property was evident in our earlier studies of the complex-temperature phase diagrams for 



the zero- field Ising model on the square, triangular, and honeycomb lattices R p4[, as well 



as the heteropolygonal 3 ■ 12 and 4 ■ 8 lattices |23 . 

Finally, we can show a close connection between the properties of the theory in the 
vicinity of the origin in the p plane by a similar argument. Small negative p, i.e. p —>■ 0~, 
is equivalent to small positive p, p —>■ + , which, in turn, corresponds to infinitely strong 
external field, h — ► 00. These results are easily seen as follows. Let us write p = p At e* 6l ' J . 
Then from the definition fi = e~ 2h , 



Y 1 

-ln/x = --(lnp M + i^ + 27rm) (2.2.8) 



h = - 

2 

As — > 0, h — > 00. In this limit, the term involving the argument 9^ has a negligible 
effect. (This actually yields a stronger result than we need here, viz., that as — > along 
any direction, not just 9^ = or tt, this direction becomes asymptotically unimportant, and 
the effect is like a physical uniform field.) This implies that in the limit as p — ► (along 
any direction), the complex-temperature phase diagram in the u plane consists only of the 
FM phase, with the AFM phase frozen out. As a consequence of the h — > —h, p —>■ 1/p 
symmetry mentioned above, this also applies (with M — > —M) to the limit \p\ — > 00 along 
any direction. 

3 Complex- Temperature Properties for Solved Cases 
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3.1 h = (/i= 1) 

In order to achieve a better understanding of the results that we shall obtain concerning 
the complex-temperature phase diagram for the interval — 1 < fi < 1 (or equivalently, the 
union {— oo</i<— 1} U {l< / u<oo}) where the model has not been solved exactly, 
it is useful to discuss this phase diagram for the two cases where it is known precisely, viz., 
fi = ±1. As part of Fig. 1, we show the complex-temperature phase diagram for H = 
(fi = 1). The phase boundaries in the u plane form a limagon ||, given by 

Re{u) = 1 + 2 3/2 coscj + 2cos2w (3.1.1) 

Im(u) = 2 3/2 sin uo + 2 sin 2oo (3.1.2) 

traced out completely for < uo < 2tt (see Fig. 1). The complex-temperature extensions 
of the physical phases with spontaneous Z 2 symmetry breaking via ferromagnetic and an- 
tiferromagnetic long range order lie, respectively, within the inner branch of the limagon 
and outside the outer branch of this limagon. The complex-temperature extension of the 
Z2-symmetric, paramagnetic (PM) phase lies between the inner and outer branches of the 
limagon. Note that the FM, AFM and (two wedges of the) PM phase are all contiguous 
at the point u = — 1 = u s . Recall that the physical critical point separating the PM and 
FM phases is u c = 3 — 2y / 2 = 0.171573..., and, in accordance with the u — > 1/u sym- 
metry noted above, the corresponding critical point separating the PM and AFM phase is 
u = l/u c = 3 + 2y / 2 = 5.82843.... Henceforth, we shall generally refer to these phases 
simply as FM, AFM, and PM, with the qualifier "complex-temperature extension" being 
understood. The corresponding phase boundaries in the z plane consist of the intersecting 
circles z = ±1 + 2 1 ^ 2 e ld , < 9 < 2n [^, |J. (In this variable, there are actually four phases: 
FM, AFM, PM, and O, where the O phase has no overlap with any physical phase ||.) 
The spontaneous magnetization is 

M(u, h = 0) = K - ; (i_ l)i/2 ( 3 - L3 ) 

in the physical FM phase 0, and this expression holds, by analytic continuation, throughout 
the complex-temperature extension of the FM phase (M = elsewhere). Thus, M vanishes 
continuously both at the physical critical point u = u c and at the point u = — 1 = u s . M 
vanishes discontinuously elsewhere along the boundary of the (complex-temperature exten- 
sion of the) FM phase. As discussed in Ref. ||, the apparent singularities at u = l/u c and 
u = 1 do not actually occur since these points lie outside the FM phase where the above 
expression applies; indeed, M vanishes identically in the vicinity of u = l/u c and u = 1. By 
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a well-known symmetry, the staggered magnetization is given by M st (y) = M(u — > y), where 
y = l/u, in the AFM phase and is zero elsewhere. In addition to its well-known divergence 
at the physical critical point, the susceptibility x diverges at u = — 1 with exponent inferred 
from series analysis to be 7^ = 3/2 |Tj| |5|. Ref. || obtained the relation 

% = 2(7-1) (3.1.4) 
which explains the value of 7^ in terms of that of the usual exponent 7 = 7/4. 

3.2 h = m/2 (/1 = -1) 

In Ref. |T]| we determined the locus of points across which the free energy is non-analytic, 
and the corresponding complex-temperature phase diagram, for the Lee- Yang solution at 
fi = — 1. These consist of the union of the unit circle and a certain line segment on the 
negative real u axis: 

{ u = e ie , 9 E [0, 2tt) } |J { l/u e < u < u e } (3.2.1) 

where 

u e = -(3-2V2) (3.2.2) 

(which is also equal to minus the value of the critical point u c in the h = model.) Since 
the nonzero value of H explicitly breaks the Z 2 symmetry, there is no PM phase. The 
FM and AFM phases occupy the interior and exterior of the unit circle, respectively. In 
contrast to the h = case, here a subset of the continuous locus of points across which / 



is non-analytic, viz., the line segment in ( [3.2.1 ), does not completely separate any phases, 



but instead protrudes into the FM phase and AFM phases, terminating in the respective 
endpoints u = u e and u = l/u e . The magnetization is 

fl +u) 1/2 

M[u, h = in/2) = f— -— (3.2.3) 

(1 - uyl^l - u/u e )(l - U e U)\ L '* 



within the FM phase (and zero elsewhere). As is clear from Q3.2.3 ), this vanishes continuously 



at u — — 1 with exponent f3 s = 1/2, and diverges at u = u e with exponent (3 e = —1/8 and 
at u = 1 with exponent (3\ = —1/4. Elsewhere on the boundary of the FM phase, i.e., the 
unit circle in the u plane, M vanishes discontinously. Note that the apparent divergence at 
the point u = l/u e does not actually occur, since this is outside of the complex-temperature 
FM phase, where the above analytic continuation is valid. From analyses of low-temperature 
series, we concluded that % h as divergent singularities (i) at u = u e with exponent 7g = 5/4, 
(ii) at u = 1, with exponent 7J = 5/2, and (hi) at u = u s — — 1, with exponent 7^ = 1 JL9J. 
(The actual values obtained from series analysis were 7g = 1.25 ±0.01, 7^ = 2.50 ±0.01, and 
7^ = 1.00 ±0.08.) 
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3.3 Partition Function Zeros for Solved Cases 
3.3.1 /i = l 

Since one of our two main methods of gaining information about the complex-temperature 
phase diagram of the model for \x 7^ ±1 is the use of complex-temperature zeros of Z, it is 
important to see how accurate this method is for the two cases where one knows the phase 
diagram exactly, viz., // = ±1. That is, we wish to ascertain how well, for finite lattices with 
various boundary conditions, the pattern of zeros resembles the locus of points across which 
/ is non-analytic in the thermodynamic limit. 

In order to compute the zeros, we calculate Z for finite lattices with specified boundary 
conditions (BC's). We have done this by means of a transfer matrix method ||25|| , and have 
used both periodic and helical boundary conditions (respectively PBC, HBC). We recall 
that in order to avoid frustration of (short- or long-range) AFM ordering, it is necessary and 
sufficient that the lengths L\ and L2 of the lattice be (i) both even for PBC and (i) one even, 
the other odd, for HBC. We have incorporated this restriction in our work. 

In Fig. 1 we show the comparison for the h = case. We have found that for similar-size 
lattices, the calculation with helical boundary conditions yields zeros which lie generally 
slightly closer to the limagon than the calculation with periodic boundary conditions. Ac- 
cordingly, we show in Fig. 1 the zeros for a 7 x 8 lattice with helical boundary conditions. 
In the vicinity of the FM-PM and PM-AFM critical points the density of zeros calculated 
on the finite lattice decreases as one approaches the real axis, in accord with the exact re- 
sult 0, |10| that in the thermodynamic limit this density g vanishes, in our notation, like 
g(u,n= 1) ~ |(1 — 11/ it c ) (1 — uu c ) I as u — > u c or it — > l/u c along the curves of zeros, where 
a = a' is the specific heat exponent, for both of these critical points. As one approaches the 
singular point u — — 1 from within the FM, AFM, or PM phases, the specific heat C has 
logarithmic divergences |J, i.e., a s = a' s = Together with a generalization of the above 
relation for g, viz., 

g ~ |1 - u/using] 1 ^ ^- as u — > u sing . (3.3.1) 

along the curve(s) where g has nonzero support, this implies that g(u, /x = 1) ~ |l+w| — > 
as«^ —1 along the curves of zeros comprising the complex-temperature phase boundaries. 
The behavior of the zeros in Fig. 1 is evidently consistent with this exact result. Indeed, one 
observes that the zeros avoid the region u = — 1. A similar tendency of the zeros not only 
to become less dense but also to deviate from the limagon can be seen in the vicinity of the 
PM-AFM critical point, where the last complex-conjugate pair of zeros occur significantly 
within the curve. (By the u — > 1/u symmetry, this is equivalent to the inverse deviation 
of the last complex-conjugate pair of zeros near the FM-PM critical point.) The fact that 



S 



there are no zeros on the positive real u axis is, of course, also implied by the property that 
for H = 0, the partition function for a finite lattice is a generalized polynomial in u with 
positive coefficients. 



3.3.2 fx = -l 

We have performed a similar comparison for [i = — 1 and show a plot of the zeros in Fig. 2, 
again for a 7 x 8 lattice with helical boundary conditions.^ One sees that all of the zeros in 
the right-hand half plane, and several other zeros, lie exactly (to within numerical accuracy) 
on the unit circle. The zero at u — 1 has multiplicity N s /2 (= 28 here), consistent with the 
fact that in the thermodynamic limit, the (reduced) free energy / contains an additive term 
(1/2) ln(l — u) (see Ref. || and eq. (2.29) of Ref. |Jl9|)0 corresponding to a delta function in 
the density of zeros, g, at u — 1. Because of this delta function singularity in g, the relation 
expressing the specific heat as an integral of g over the region where it has nonzero support 



receives an additive contribution. In turn, this has the effect of modifying the relation ( |3.3.1| ) 
connecting g with the critical exponent in the specific heat at a singular point u S i ng ^ Hence, 
the fact, shown in Ref. |fL9j , that the specific heat C has a finite non-analyticity (a[ = 0) at 



u = 1 as this point is approached from the FM or AFM phases does not conflict with the 
non- vanishing of g at u — 1. On the negative real axis, there are zeros lying on (or near) 
the line segment in (|3.2.1|) . As one can see from Fig. 2, near the inner (and by the u — > 1/u 
symmetry, also the outer) endpoint of the line segment, the density of zeros approaches a 
nonzero constant. This finite-lattice feature is in good agreement with the infinite-lattice 
relation ( |3.3.1|) between g and the specific heat singularity, together with the exact exponent 



which we found earlier | 19|| , a' e — 1, at u — u e and u = l/u e . Proceeding next to the 
intersection point at u = — 1, we observe that our exact result that the specific heat has a 
finite non-analyticity (a' s = 0) at u = — 1 as this point is approached from either the FM or 
AFM phases ]19| implies that in the thermodynamic limit, g(u, n = —1) ~ — » as one 



approaches this point along the unit circle or the negative real axis. This is consistent with 
the results in Fig. 2. (It happens that for the 7x8 lattice with HBC there is a zero precisely 
at u = — 1, but this does not occur for other lattice sizes and is not a general feature.) Finally, 
we note a contrast with the h = case: there, Z was a generalized polynomial in u with 



2 To save space, we truncate the plot on the left at Re(u) = —2; there are four zeros lying on the negative 
real axis, and a pair lying near this axis, to the left of this point, which are not shown. However, the plot as 
shown contains all the information, since these six zeros are just the inverses, by the u — * 1/u symmetry, of 
the six which are shown, lying on or near the negative real axis within the unit circle. 

3 As discussed in Ref. |l£j, in addition to this additive term there is also a subdominant contribution to 
the singularity in / at u = 1 arising from the vanishing of the argument of the logarithm in eq. (2.29) of 
that paper. 



9 



positive coefficients, whereas here the coefficients are not all positive, so that zeros can, and 
indeed do, occur on the real u axis. One quantitative measure of how well the zeros on the 
finite lattices approach the locus of zeros on the infinite lattice is provided by how close the 
innermost zero on the negative Re(u) axis is to the infinite-lattice result, u e = — (3 — 2y / 2)- 
We find that the position of this zero is (i) u = —0.1822 on an 8 x 8 lattice with PBC, and 
(ii) u = —0.1831 on a 7 x 8 lattice with HBC. These points are respectively 6.2 % and 6.5 
% farther away from the origin than u e . 

4 Complex- Temperature Properties < h < oo (0 < 

/!< 1) 

4.1 Information from Zeros of Z 

Having these comparisons with exactly solvable cases as background, we now proceed to 
investigate the complex-temperature phase diagram and associated singularities for the Ising 
model in a nonzero magnetic field, i.e., < H < oo, or equivalently, < /i < 1. Since 
any nonzero magnetic field explicitly breaks the Z 2 symmetry of the Hamiltonian, the Z 2 - 
symmetric, PM phase immediately disappears. We recall from section ( |2.2| ) that H can be 
taken positive without loss of generality, and that the results for the above interval in fi 
also describe the behavior for 1 < \i < oo. Of the four symmetries in section ( j2.2| ), the 
u — > 1 ju symmetry of the phase diagram (with corresponding interchange of FM and AFM 
phases) ceases to hold if fi ^ ±1. One may inquire how the complex-temperature phase 
boundary changes as H is increased from zero. One knows rigorously that for h ^ 0, it 
is possible to analytically continue from u = out along the positive real u axis through 
the point u c (which is no longer a critical point), all the way to the vicinity of u — 1 
(K = 0). This means that the region which was previously, for H = 0, the PM phase 
and its complex-temperature extension, immediately becomes part of the FM phase and its 
complex-temperature extension, respectively. This also means that for any nonzero h, the 
inner branch of the limagon which formed the phase boundary between the FM and PM 
phases for H = breaks apart at u = u c into two separate arcs, leaving an opening through 
which the analytic continuation mentioned above can be performed. 

One also knows that if J < 0, then for a fixed value of H, at sufficiently low T, the 
system will still exhibit long-range AFM ordering. Equivalently, for fixed h and hence fi, 
and sufficiently large u > 1 (i.e. sufficiently large negative K), there will be AFM ordering. 
Thus, for any finite h, there is still a non-analytic boundary separating the physical and 
CT-extension of the FM phase from the physical and CT extension of the AFM phase. Since 
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the nonzero H biases the system toward FM ordering and hence against AFM ordering, it 
follows that the singular point which separates the physical FM and AFM phases, u FM AFM , 
increases as H increases (starting from the value liin.H_ > o u PM AFM = l/u c ). As u increases 
through the value u FM AFM , the system undergoes a first-order transition at which M jumps 
discontinuously to zero. One expects that generic points on the border between the (complex- 
temperature) FM and AFM phases are also associated with a first -order phase transition. 
Given the biasing effect of nonzero H, one expects the complex extension of the FM phase 
to increase as H increases. Because of the symmetry under h — > —h noted above, it follows 
that the rightward shift in u FMAFM as a function of H is independent of the sign of H. 
An additional general statement is that for all H (zero or not), the (complex-temperature 
extensions of the) FM and AFM phases must be completely separated by a phase boundary. 

In Fig. 3 we show a combined plot of the CT zeros of Z, calculated for the same lattice 
(7x8 with helical boundary conditions) as before, for a sequence of increasing values of 
h, and in Fig. 4 we show plots of these zeros for certain fixed values of h. There are 
several interesting features which we can discern from these calculations. The results are 
in accord with the expectations based on rigorous arguments noted above. As h increases 
above zero (// decreases below 1), the right hand end of the inner branch of the limagon which 
formerly crossed the real u axis at u c immediately breaks at u = u c , rendering it possible to 
analytically continue from u = along the positive real axis outward past u — all the way 
to the vicinity of l/u c . This breaking yields two complex-conjugate arcs, which retract from 
the real u axis as h increases. Scaling arguments imply that in the neighborhood of h — 0, 
the endpoint of the arc in the upper half plane should move in the direction [^] 

«- = ^ ("-I) 

For the 2D Ising model, inserting the exactly known critical exponents (3 = 1/8 and 5 = 15, 
it follows that 

0ae,2D = ^ = 4&° (4.1.2) 

One can see roughly, and we have confirmed accurately that this is in agreement with the 
motion of the zero nearest to u c for small h. The direction of motion away from the origin 
is different for the zeros on the inner arcs. Related to this, as h increases, the upper arc 
becomes oriented in a more northwest to southeast, rather than west to east, direction. 

A notable feature that one sees in each of the graphs in Fig. 4 is that there is a sharp 
increase in the density of zeros as one approaches the endpoints of these arcs. This suggests 
the possibility that in the thermodynamic limit the density g might actually diverge as one 
approaches these endpoints along the arcs. Below, we shall find support for this conclusion 
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from our analysis of low-temperature series, which yields a value for the specific heat exponent 
a' ae > 1 at these endpoints; together with the relation ( |3.3.1|) , such a value of a' ae implies 
that g diverges at the arc endpoints. Note also, of course, that a' ae > 1 also means that the 
internal energy diverges at these endpoints. In this context, we recall that in our study of 
the complex-temperature singularities of the Ising model on the honeycomb lattice |24| , we 
found the exact result that at the point z = —1 = Z£, as approached from either the FM or 
AFM phases, the specific heat diverges with an exponent a' e = 2. 

In accord with the argument given above, as h increases, u FMAFM increases, and, more 
generally, one sees that in the right-hand half plane, the zeros on what was the outer branch 
of the limagon for h = move outward. The directions of motion of these zeros are approxi- 
mately radially outward and do not change significantly as h is increased. Of course, the fact 
that the zeros on what were, for h = 0, the inner and outer branches of the limagon both 
move outward is a striking manifestation of the loss of the u — > 1/u symmetry which held 
for h = 0, since the latter symmetry implied that these zeros were inverses of each other, 
(so that if the inner branch moved out, the outer branch would have had to move in). From 
the comparison of the zeros calculated on finite lattices with the exact phase boundary for 
H = 0, we saw that the pair of zeros farthest to the right lies within the outer branch of 
the limagon. One thus anticipates that this may also be true for nonzero H. Numerical esti- 
mates of K b = J I (k B T b (H)), expressed in terms of the u variable, are u FM AFM = 5.90909(9), 
6.34036(9), and 7.99333(10) for h = 0.2, 0.5, and 1.0, respectively, where the number in 
brackets is the uncertainty in the last digit |§. If we draw a curve through the zeros on the 
right-hand side of the phase diagram, and then compare the point at which they cross the 
real axis with these values, we see that, indeed, the finite lattice calculations yield a slight 
underestimate of this point for nonzero H, as they did for H = 0. 

The motion of the zeros in the left-hand half-plane is more complicated. One sees that, 
while the zeros do move generally outward from the origin as H increases, the directions of 
motion of a number of zeros undergo significant (and non-monotonic) changes. We note that, 
owing to the loss of the inversion symmetry for fi ^ ±1, the self-inverse point u — — 1 loses 
the special significance which it had for these two values, for which it was an intersection 
point on the respective phase boundaries. 

From these calculations, we may draw some further inferences about the phase diagram 
in the thermodynamic limit. As h increases, the boundary between the FM and AFM phases 
on the right-hand side of the phase diagram moves outward from the outer branch of the 
limagon. On the left-hand side, the intersection point at u = — 1 is replaced by an outer 
boundary to the left of this point, which shows some concavity but moves farther to the left 
as h increases. Recall that for h = 0, the inner and outer branches of the limagon cross in a 



12 



perpendicular manner at this point || (see Fig. 1). Our calculations suggest that it is likely 
that as h increases from 0, the curves which, for h = 0, were the inner and outer branches 
of the limagon, no longer intersect at u = —1; instead, the upper and lower parts of the 
former outer branch join smoothly to each other on the left, and cross the real u axis in a 
vertical manner rather than as perpendicular NW-SE and SW-NE curves. Our results also 
suggest that the two complex-conjugate arcs which, for h = 0, comprised the connected inner 
branch of the limagon, no longer cross the negative real u axis at u = — 1 but instead retract 
from this axis. It is possible that they join the outer phase boundary (at complex-conjugate 
points). As h increases more, our calculated zeros indicate that these complex-conjugate 
arcs move progressively farther from the real axis. As h — > oo, the boundary separarating 
the FM and AFM phases moves outward to complex infinity in the u plane, so that in this 
limit, the entire plane is occupied by the FM phase. We shall next obtain further information 
about the endpoints of these arcs from a series analysis. 



4.2 Results from Series Analysis 

In order to investigate the complex-temperature singularities for — 1 < /i < 1, we shall make 
use of the low-temperature, high-field series expansion for the free energy of the Ising model 
on the square lattice |E?J-[23|. The partition function Z can be written as Z = e Na ^ 2K+h ' Z r , 



and equivalently, the reduced free energy can be written as / = 2K + h + f r , where 
f r = lini7v s ^oo N' 1 In Z r . In Ref. [E3], Baxter and Enting calculated the low-temperature 



expansion of the quantity k = e^ r = limA^oo Z\l Ns to 0(u 23 ) (with the coefficients of the 
powers of u being exact polynomials in n). This expansion has the form 

oo 

«: = l + EE«n, m M> m (4.2.1) 

n=2 m 



where j < m < j 2 for n = 2j and j < m < j(j — 1) for n = 2j — 1. From this series, 
we calculate, for each value of fi, resultant series for the specific heat C, magnetization M, 
and susceptibility x- We then analyze these using dlog Pade and differential approximants 
(abbreviated PA and DA; for a recent review, see [^9|). Our notation for differential approx- 
imants follows that in Ref. |29] and our implementation is the same as that in Ref. [ 19[ ; in 
particular, we use unbiased, first-order differential approximants. 

We find from both the PA and DA analyis that for < fi < 1, the series indicate a 
complex-conjugate pair of singularities at positions which are in very good agreement with 
the locations of the innermost zeros on the inner arcs. We denote these arc endpoints, u ae 
and u* ae . Table 1 shows the comparison of the values of the arc endpoint u ae as obtained from 
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h 




U ae ,Z 


Uae,ser. 




ftae 




0.8 


0.202 


0.313 + 0.663i 


0.3203(6) + 0.6503(6)i 


1.175(7) 


-0.18(5) 


1.195(10) 


1.0 


0.135 


0.320 + 0.857i 


0.3265(5) + 0.8426(6)i 


1.177(8) 


-0.18(5) 


1.20(2) 


1.2 


0.0907 


0.324 + 1.086i 


0.3306(5) + 1.0685(6)i 


1.18(2) 


-0.18(5) 


1.19(2) 



Table 1: Values of the arc endpoint, u ae , for several values of h, from (i) calculation of zeros 
of Z on finite lattices, yielding u a£) z, and (ii) analyses of series for C, M and x, yielding 
u ae ,ser.- Table also gives values of the exponents a' ae , /3 ae , and •y' ae at u ae as determined 
from analysis of low-temperature series using dlog Pade and differential approximants. The 
values of u ae , S er. are quoted to an accuracy in accord with the agreement between the different 
series and methods of analysis. The values for u ae> z are quoted to an accuracy reflecting the 
differences in values with different boundary conditions. The numbers in brackets denote 
the uncertainties in the last digit of a given entry. 



the complex-temperature zeros on finite lattices and from the analysis of the low-temperature 
series. The fractional differences in the positions are about 2 %. 

We have made a preliminary study of the singularities at these arc endpoints. Since the 
arc endpoints arise abruptly with the breaking of the inner branch of the limagon for any h, 
regardless of how small, one anticipates a similarly abrupt change in exponents. Indeed, in 



our previous studies of complex-temperature singularities (e.g. Refs. |24], p3|), using exact 
results, we found that M always diverges at endpoints of arcs protruding into the FM phase. 
Hence, we expect (and will verify) that M diverges at the arc endpoints in the present case, 
so that f3 jumps discontinuously from the value (3 = 1/8 for h = to negative value(s) for 
nonzero h. We assume the usual leading singular forms 

C~A c (l-u/u ae )- a '<- (4.2.2) 

M ~ A M {1 - u/uaef™ (4.2.3) 

X~A x (l-u/u ae )-^ (4.2.4) 

(We do not consider confluent singularities here; at least for H = 0, they are unimportant 
for the 2D Ising model, a fact used also in our previous study of complex-temperature 
singularities in this case 0.). We analyze the singularities for three values of h which are 
chosen to be sufficiently large that the arc endpoints are reasonably well separated from each 
other, viz., h = 0.8, 1.0, and 1.2. We obtain the results shown in Table 1. In the case of a' ae 
and 7^ e , the dlog Pade and differential approximants yielded values in very good agreement 
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with each other. In the case of f3 ae , the PA's gave somewhat larger values than the DA's; 
we list values which are more heavily weighted by the (presumably more accurate) DA's, 
with estimated errors which are increased to reflect the PA results. As is evident from Table 
1, for the values of h which we have considered, these exponents are consistent with being 
independent of the value of h in this range. The values of the exponents are also consistent 
with the equality a' ae + 2f3 a e + l' ae = 2. The fact that j3 ae is negative means that M diverges at 
the arc endpoints; this is the same behavior that we found earlier in our studies of complex- 
temperature singularities on various lattices, using the respective exact expressions for the 
magnetization (e.g. see the summary in Table 2 of Ref. f2"3fl). Our results also indicate that 
a' ae > 1. As noted above, in conjunction with the relation ( p. 3.1 ), this implies that in the 



thermodynamic limit the density g of zeros of Z diverges as one moves along the arcs and 
approaches their endpoints. The behavior of the zeros calculated on finite lattices already 
suggests this infinite-lattice result. 

5 Phase Diagram for — 1 < fi < 
5.1 Zeros of Z 

This range of fi is not physical, in contrast to the range < fi < 1 discussed above. 
Nevertheless, it is of interest because by determining how the phase diagram changes as 
\i decreases through this range, we can complete the connection between the exact solutions 
at fi = 1 and —1. 

For the discussion of the complex-temperature zeros of Z, just as in the previous section, 
we started with the exactly known case /i = 1 and moved away from it by decreasing fi, so 
also it is convenient here to start with the exactly known phase diagram for fi — — 1 and 
move away from it by decreasing the magnitude of fi toward 0. In Fig. 5 we show a plot 
of the paths of individual zeros for a sequence of values of h, and in Fig. 6 we show plots 
of zeros for fixed values of h. These results have a number of salient features. First, as 
fi moves from —1 toward 0~, the zeros which were on the unit circle move approximately 
radially outward to larger values of \u\. They maintain a roughly circular form. The zeros 
which had been on or close to the negative real axis and had corresponded to what, in the 
thermodynamic limit, was the line segment in ( |3.2.1| ) extending between l/u e and u e for 



fi = —1, move gradually outward, away from the origin. Meanwhile, for any value of fi 
greater than —1, there abruptly appears a new set of zeros on the positive real axis. Our 
calculations of zeros on finite lattices suggest the inference that in the thermodynamic limit, 
these form a finite line segment of singular points starting initially at u = 1 for fi = — 1 + e, 
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e — > (and extending rightward). As fi increases from —1, the right-hand phase boundary 
between the FM and AFM phases moves outward sufficiently quickly that, as one can see 
in Fig. 6, it lies to the right of the new line segment of zeros on the positive real u axis. 
We shall denote the respective endpoints, in the thermodynamic limit, of the left-hand and 
right-hand line segments of zeros as uih e and u r h e - A striking feature of the plots is that 
the density of zeros increases sharply as one moves rightward toward the inner end of the 
line segment of zeros on the negative real u axis, and also as one moves leftward toward 
the inner end of the line segment of zeros on the positive real axis. This is quite similar 
to what we observed for the density of zeros near the arc endpoints in the FM phase for 
< fi < 1. Again, it suggests the inference that in the infinite-lattice limit, the density 
of zeros g diverges at these endpoints uih e and u r h e - Our series analysis below support this 
inference. The results in Fig. 6 also show that as fi moves toward through the interval 
— 1 < fi < 0, the endpoints of the line segments, uih e and u r h e move slowly away from the 
origin. As fi — > _ , the boundary separating the FM and AFM phases moves outward to 
complex infinity, so that in this limit, only the FM phase remains. This can be understood 
on general grounds, as we have discussed above. 



5.2 Series Analysis 

From our analysis of the low-temperature series for C, M, and x, we have obtained very 
good agreement with the positions of the inner endpoints of the left and right line segments, 
u lhe and u rhe . Table 2 displays this comparison. Of course, one cannot use the series to 
locate the outer endpoints of the left and right line segments since these lie outside the FM 
phase where these series are applicable. 
Since 

lim u ihe = u e (5.2.1) 



and since the exponents are known for \x = —1 as |B| a' e = 1, (3 e = —1/8 (exact), and 
7g = 5/4 (inferred from series analysis), an interesting question is whether, for values of \i 
in the interval — 1 < /i < 0, the singularities in C, M, and x & t the inner endpoint of the 
left line segment are described by exponents which are (i) independent of \x in this interval, 
and (ii) equal to their values at fi = — 1. In Table 3 we list our results for \i = —0.5 
and —0.2 (the results for /i = —0.1 were less precise). The values for a' lhe and 7 ; ' he from 
the dlog Pade and differential approximants agreed very well with each other. The value 
of fyke was more difficult to extract from the series analysis (in particular, the dlog Pade 
approximants yield somewhat larger values than the differential approximants); using only 
the (presumably more accurate) DA's, we obtain the value given in Table 3. To within the 
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h r 


n 


U lhe,Z 


lne,ser. 


U rhe,Z 


u h 

rne.ser. 


0.347 


-0.5 


-0.2976 


-0.2870(5) 


1.047 


1.053(15) 


0.805 


-0.2 


-0.5467 


-0.5314(8) 


1.246 


1.245(10) 


1.151 


-0.1 


-0.8471 


-0.8267(8) 


1.529 


1.522(10) 



Table 2: Values of the inner endpoints of the left and right line segments for negative fx, 
with h = — (1/2) In /i = h r + m/2. Entries in the columns denoted Uih e ,z and u rheZ are 
the innermost zeros on the left- and right-hand line segments, calculated on a 7 x 8 lattice 
with helical boundary conditions. Entries in the columns denoted Uih e ,ser. and u r h e ,ser. are 
the positions of the inner endpoints of the left and right line segments as obtained from the 
analyses of low-temperature series for C, M, and % Accuracies listed are chosen as in Table 
1. 



n 


a 'lhe 


Plhe 


llhe 


a 'rhe 


ftrhe 


^rhe 


-0.5 


1.20(5) 


-0.2(1) 


1.20(4) 


1.65(12) 


-0.35(20) 


1.00(15) 


-0.2 


1.18(5) 


-0.2(1) 


1.19(5) 


1.4(1) 


-0.2(1) 


1.15(20) 



Table 3: Values of the singular exponents a' lhe , flihe, lihei a 'rhei Prhe, and ^' rhe at the inner 
endpoints of the left- and right-hand line segments, u\h e and u r h e , respectively, as determined 
from analyses of low-temperature series for C, M, and x using dlog Pade and differential 
approximants. 

estimated uncertainties, these results are consistent with being independent of fi for the two 
values of fx considered, and with satisfying the equality a' lhe + 2/5^ e + 7^ e = 2. Regarding 
question (ii), 7^ e is slightly, but not decisively, lower than the inferred value at \x = — 1, 
7g = 5/4. The exponent f3ih e is also consistent with the fx = — 1 value, f3 e = —1/8, but the 
large uncertainty in its determination precludes a sensitive test of this equality. The values 
of a' lhe are about 4<r above unity. Assuming the series analysis gives a reliable determination 
of this exponent, this means that it is not equal to the fx = — 1 value, a' e — 1. Furthermore, 
it means that the density of zeros does diverge at w^ e , a conclusion already hinted at by the 
behavior of the zeros calculated on finite lattices. 

We have also extracted the critical exponents at the inner endpoint of the right-hand 
line segment u r h e from series analysis. Here there are two specific questions to investigate. 
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First, are these exponents consistent with being independent of fi in the range —1 < fi < 0? 
Second, recalling that, in the limit fi — > — 1, the right line segment decreases in length and 
eventually disappears at u — 1 (as inferred from our finite-lattice calculations, in agreement 
with our series analysis), one would like to see whether the exponents on the inner endpoint of 
this right-hand line segment, u r h e , might be equal to those which we found for the singularity 
at u = 1 for fi = — 1, viz., a[ = (finite logarithmic non-analyticity in C), (3\ = — 1/4 (both 
exact), and j[ = 5/2 (inferred from series analysis) |19[ . Our results are shown in Table 3. 
Just as was true of the inner endpoint of the left line segment, it was more difficult to get an 
accurate value of p r h e , and again we have quoted a value weighted more by the differential 
approximants, with a commensurately large uncertainty assigned. These values of critical 
exponents at the inner endpoint of the right-hand line segment are crudely consistent, to 
within large uncertainties, with being independent of fi over the indicated range of fi. They 
are also in accord with the hypothesis that a' rhe + 2j3 rhe + j' rhe = 2. As fi — > — 1, we find that 
it is more difficult to get accurate values for the singular exponents. This is not surprising, 
since in this limit the right line segment disappears. Regarding the second question, although 
the uncertainties in the critical exponents are large, our results indicate that, for the range 
of fi considered, the exponents a' rhe and j' rhe are not equal to the respective exponents at 
u = 1, fi = — 1. We cannot draw a firm conclusion concerning whether over the same range 
of fi, P r he is equal or unequal to f3\. The value(s) of a' rhe is about 4-5 a above 1, which 
supports the conclusion that the exact value of this exponent is, indeed, > 1, and hence that 
the density g of zeros diverges at u r h e , as it does at u\h e - This is again consistent with the 
observed behavior of the zeros calculated on finite lattices. 



6 Comments on Other Lattices 
6.1 Honeycomb Lattice 

The complex-temperature phase diagram of the zero- field Ising model on a honeycomb lattice 
was determined in Ref. from exact results. In the complex z plane the boundaries consist 
of the union of an arc of the unit circle z = e* 9 , for 71/ 3 < 9 < 5n/3 with a closed curve 
lying in the region Re(z) > which intersects the circle at z = ±2 and the positive real 
z axis at z c and l/z c , where z c = 2 — \^3. This phase diagram consists of three complex- 
temperature phases: (i) FM, surrounding the origin, (ii) PM around z = 1, and (iii) AFM, 
the outermost phase, extending to complex infinity. The complex-temperature singularities 
were determined exactly for the specific heat and magnetization, and were studied using 
low-temperature series expansions for the susceptibility |24}| . 
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Here we determine the complex-temperature phase diagram for the honeycomb lattice 
at jj, = — 1. One first uses the identity e 1 ^^ 2 ^" = ia n to relate the model to one with no 
magnetic field and the partition function 

Z = i N °Y<([[vn)e^<™'> r7nK ™' CTn ' (6.1.1) 

a n n 

(where we have allowed the possibility that the Kj along the three different lattice vectors are 
different). One can next associate dimers with pairs of a's at neighboring sites and then re- 
exponentiate these to express the partition function in terms of a model with the associated 



coupling Kj shifted by in/2 [17], 18 1. A complete dimer covering is provided by placing 
these dimers along one of the three lattice vectors on the honeycomb lattice. The situation 
is particularly simple, for the following reason: given that the coordinaton number is odd 
(q = 3), one can multiply by 1 = a 2 at each vertex and thereby place the dimers on every 
bond of the lattice, carry out the re-exponentiation, and thus map the model to one with 
the couplings along each of the three lattice vectors shifted as Kj — > Kj + in/2, j = 1, 2, 3. 
Since the shift is the same for each of these lattice directions, one can immediately specialize 
to the case of isotropic couplings. As a result, the reduced free energy / for h = in/2 
is given, up to an additive term in/ 4, simply by that for h = with the replacement 
K — > K + in/2, or equivalently, z — » — z or v — > 1/v. This is also true of C, M, M st , 
and x- Hence, we can directly take over our results for the zero-field case, and we find 
the complex-temperature phase diagram shown in Fig. 7. This consists of the FM, and 
AFM phases together with a phase in which M = M st = identically, whence the label 
ZM (standing for "zero uniform and staggered magnetization"). The magnetization satisfies 
M(z, h = in/2) = M(—z, h = 0), and similarly for the staggered magnetization. The points 
in the ZM phase are in 1-1 correspondence, under z — > —z, with the points in the PM 
phase of the zero-field model. The reason that we do not label it as PM is that the nonzero 
magnetic field explicitly breaks the Z 2 symmetry, so that there is no PM phase in the strict 



sense. However, as discussed in Ref. []T9| , this non-invariance is manifested only in the 
non-invariance of the constant in /4 in /, and hence does not affect derivatives such as C, 
M, and \. We note the contrast with the square lattice, where for // = —1 (and indeed, also 
for all //^ 1), there is no phase in which M and M st both vanish; instead of becoming the 
ZM phase, the PM phase disappears completely. Our determination of complex-temperature 



singularities in x m Ref- H3 can also be taken over directly to the present case, with just 



the change z — * —z. The phase diagram in v shows the interesting feature that the FM 
and AFM phases are separated by a line of singular points on the imaginary v axis. This is 
similar to what we found for the square lattice (see Fig. 2(b) of Ref. ||19|| ) (A difference is 



that in the case of the square lattice, the FM and AFM phases extend to complex infinity in 
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the right and left parts of the v plane, and hence the line separating them is the complete 
imaginary v axis, whereas for the honeycomb lattice, the FM and AFM phases are bounded 
by an oval-like curve (Fig. 7(b)) outside of which is the ZM phase, and the line segment 
separating them is finite, extending from — i\f?> to i^/3.) 

As h increases from 0, the rigorous argument that one can analytically continue outward 
from z = past the former FM-PM critical point at z c , to the vicinity of the former PM- 
AFM critical point at l/z c implies, as in the case of the square lattice, that the right-hand 
boundary of the FM phase breaks open at z = z c , forming two complex-conjugate arcs which 
retract away from the real axis. The former PM phase becomes part of the enlarged FM 
phase, again as in the square lattice case. The FM phase is completely separated from the 
AFM phase, and the former second-order critical point at l/z c which separated the PM and 
AFM phases for h = now becomes a first-order phase boundary between the FM and AFM 
phases. The general arguments given above imply that as fi — > + (and also for fi — ► CP), 
the FM phase expands to fill the entire z plane. Finally, for the special value \l — — 1, the 
model again exhibits a phase with M = M st = 0, the ZM phase. 

It may be noted that the simple relation under z — > — z between the free energy and 
other thermodynamic functions which holds for this lattice also holds for heteropolygonal 
lattices with odd coordination number. Hence the complex-temperature phase diagrams for 
the 3 • 12 2 and 4 ■ 8 2 lattices with h = m jl are simply given by the diagrams which we worked 



out for the zero- field case [23 with the replacement z — > —z (and PM — > ZM). 



6.2 Triangular Lattice 

We recall that for the zero-field Ising model on the triangular lattice, the locus of complex- 
temperature points where the free energy is non-analytic consists of the union of the circle 

„ = -i + | e * (6.2.1) 

(with < 6 < 2tt) and the semi-infinite line segment 

1 

-oo<u<-- (6.2.2) 
o 

The FM phase is located within the circle, with the PM phase lying outside and extending 
to complex infinity. 

By methods similar to those we have discussed, we calculate the complex-temperature 
phase diagram for h = iir/2 and show it in Fig. 8. In the u plane, the continuous locus of 
points where the free energy is non-analytic consists of the union of the circular arc 

u = -{-l + e ie ), e ce <\6\<7i (6.2.3) 
2 
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where 



corresponding to the endpoints u, 



e, 



ce 



ce 



arctan^-^j ~ 38.9' 

e t9 ce anc [ wnere 



(6.2.4) 



u, 



ce 



9 



(6.2.5) 



and the semi-infinite line segment 



1 

— oo < u < — - 



(6.2.6) 



Just as there was no AFM phase for h = 0, so also there is none for h — m/2. As h increases 
from 0, the PM phase abruptly disappears. By the same rigorous arguments as were discussed 
before, the right-hand boundary of the FM phase breaks open at u = u c = 1/3 and the circle 
changes into a circular arc, the complex-conjugate endpoints of which retract away from the 
positive real axis. All points in the u (equivalently, z or v) plane (not lying on the above 
continuous locus of points where / is non-analytic) are analytically connected to each other 
and all lie in the FM phase. The corresponding locus of points in the v plane consists of the 
union of a circular arc v = e'^ traced out for ce < < ir, where ce = arctan(2 3 / 2 /3) = 43.3° 
with two complex-conjugate arcs intersecting this circle at v — ±i. Comparing the h = 
and h = in/2 complex-temperature phase diagrams in the u plane, one sees that the former 
circle of radius 2/3 centered at u = —1/3 is replaced by an arc of a circle of radius 1/2 
centered at u = —1/2, and the semi-infinite line segment extending leftward from u = —1/3 
is replaced by a similar line segment starting at —1/2. 



In this paper we have determined the general features of the complex-temperature phase 
diagram for the square-lattice Ising model in a nonzero external magnetic field < h < oo 
(0 < n < 1). Our methods included calculations of complex-temperature zeros of the 
partition function for finite lattices and analysis of low-temperature series. We have also 
performed a similar analysis for — 1 < /i < 0. The results enable one to exhibit a continuous 
connection, via the variation of /i through real values from 1 to —1, of the two known 
exact solutions of this model, viz., those of Onsager and of Lee and Yang. We have studied 
the exponents in the specific heat, magnetization, and susceptibility at certain complex- 
temperature singularities which are present in the interval — 1 < /i < 1. Extending our 
earlier work, we have also determined the exact complex-temperature phase diagrams of the 



7 Conclusions 
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model at [i — — 1 on the honeycomb and triangular lattices and have discussed the relation 
between these and the zero-field phase diagrams. Our results give further information about 
the still-intriguing issue of the Ising model in a magnetic field. 

This research was supported in part by the NSF grant PHY-93-09888. 
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Figure Captions 



Fig. 1. Complex-temperature zeros of Z, calculated for h = (// = 1) on a 7 x 8 lattice 
with helical boundary conditions, as compared with the exact result, in the u = e~ 4K plane. 

Fig. 2. Complex-temperature zeros of Z, calculated for h = in/2 (fi = — 1) on a 7 x 8 
lattice with helical boundary conditions, as compared with the exact result, in the u = e~ AK 
plane. 

Fig. 3. Zeros of Z in the u plane, for h varying from to 1.25 in increments of 0.05. 
Lattice size and boundary conditions are as in Fig. 1. 

Fig. 4. Zeros of Z in the u plane, for h= (a) 0.25 (b) 0.5 (c) 0.75 (d) 1.0. Lattice size 
and boundary conditions are as in Fig. 1. 

Fig. 5. Zeros of Z in the u plane for — 1 < /j, < 0. Denoting h = in/2 + h r , h r varies 
from 0.25 to 1.25 in increments of 0.05. Lattice size and boundary conditions are as in Fig. 
1. 

Fig. 6. Zeros of Z in the u plane for /i— (a) —0.5 (b) —0.2 (c) —0.1. Lattice size and 
boundary conditions are as in Fig. 1. 

Fig. 7. Complex-temperature phase diagram for the Ising model on the honeycomb 
lattice with \x = — 1 in the variable (a) z (b) v. 

Fig. 8. Complex-temperature phase diagram for the Ising model on the triangular lattice 
with /i — — 1 in the variable u. 
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